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Abstract. One of the important aspects in recent trends in complex analysis has 
been the increasing degree of cross-fertilization between the latter and mathematical 
physics with great benefits to both subjects. Contour dynamics in the complex plane 
turned to be a meeting point for complex analysts, specialists in stochastic processes, 
and mathematical physicists. This was stimulated, first of all, by recent progress 
in understanding structures in the classical and stochastic Lowner evolutions, and 
in the Laplacian growth. The Virasoro algebra provides a basic algebraic object 
in conformal field theory (CFT) so it was not surprising that it turned to play an 
important role of a structural skeleton for contour dynamics. The present paper is 
a survey of recent progress in the study of the CFT viewpoint on contour dynamics, 
in particular, we show how the Witt and Virasoro algebras are related with the 
stochastic Lowner and classical Lowner-Kufarev equations. 



1. Introduction 

Conformal field theory (CFT) in two dimensions has deep intrinsic connection 
to representation of infinite dimensional algebras. The Virasoro algebra is a vertex 
algebra, which appeared in early 1970's physics papers (see, e.g. [48j) on string 
theory. Earlier in 1968 it was introduced by Gelfand and Fuchs |15] as a unique 
(up to isomorphisms) central extension of the algebra of vector fields on the unit 
circle. Later in 1980's, it became clear that the Virasoro algebra turned out to be 
a universal symmetry algebra in two-dimensional CFT. The infinitesimal conformal 
transformations in the classical setup lead to an infinite dimensional algebra, called 
the Witt algebra. Turning to quantum field theories, the conformal anomaly, or Weyl 
anomaly, leads to the appearance of a nontrivial central charge. So the Witt algebra 
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is modified by central extension to the Virasoro algebra. Infinite-dimensional algebras 
have also been used recently in the theory of exactly solvable models. For example, 
the Virasoro algebra plays a central role in the study of integrable systems, such 
as those associated to the KdV and other soliton hierarchies. The Virasoro algebra 
is intrinsically related to the KdV canonical structure where the Virasoro brackets 
become just the Magri brackets for the Miura transformations of elements of the 
phase space of the KdV hierarchy (see, e.g., [101 116j). The lattice Virasoro algebra 
appears in the study of the Toda field theory and Toda integrable systems |11L I21j . 

On the other hand, contour dynamics is a classical subject in complex analysis. 
One of the typical dynamics started from classical Hele-Shaw experiments in 1897. 
This leads to a sample free boundary problem, known also as the Laplacian growth in 
two dimensions (see |20LI47j ). Another one is the Lowner evolution, see [271.140] . Last 
decade, the progress in the study of Laplacian growth has resulted in its integrable 
structure in terms of Richardson's moments, namely these moments form a disper- 
sionless Toda hierarchy [34] . Richardson's moments are conserved quantities of the 
Laplacian growth. Recently, it was proved |33j that the Virasoro generators in their 
co-vector representation are conserved quantities of the Lowner- Kuf are v evolution. 
Friedrich and Werner [13j, and independently Bauer and Bernard [6j found relations 
between SLE (stochastic or Schramm-Lowner evolution) and the highest weight rep- 
resentation of the Virasoro algebra. All these results encouraged us to conclude that 
the Virasoro algebra is a common structural basis for these and possibly other types 
of contour dynamics and we decided to present our overview of the development in 
this direction. 

The survey is designed for a comprehensive reading of sufficiently wide mathe- 
matical audience, in particular graduate students. Therefore, we decided to include 
several parts well-known for specialists in mathematical physics but less known for 
analysts. 

We start to show how the Virasoro algebra appears in CFT and in integrable 
models. Then we proceed with the classical Lowner-Kufarev equation where the 
central extension does not appear and we work only with the Witt algebra. The 
central extension appears for the stochastic version of the Lowner equation. Finally, 
we briefly describe connections between SLE and CFT. 

This survey is partially based on a lecture given by the second author at the 
II Winter School in Complex Analysis and Operator Theory, held in Seville (Spain), 
February 5-9, 2008. He is thankful to the organizers for their hospitality. The authors 
want to express their gratitude to Helene Airault, Roland Friedrich, Paul Malliavin, 
Yurii Neretin, and Dmitri Prokhorov for many fruitful discussions on the subject 
during last years. 
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2. CFT and Virasoro algebra 

There is a vast amount of standard references on basics of CFT, and classical and 
quantum strings, see e.g., |19L 137] . but since the topic is presented in a lecture form 
for mathematics students, we take a risk to look briefly through a simplest example 
of a closed bosonic string with periodic boundary conditions in order to introduce the 
Virasoro algebra. 

2.1. Classical bosonic string. We start with some basic definitions. A p-brane 
is spatial evolution of a p-dimensional object, which is reduced to a point particle if 
p = 0, a. string if p = 1, a membrane if p = 2, etc., in some D- dimensional ambient 
spacetime Tt. The result of such evolution is called the worldhne (p = 0), worldsheet 
(p = 1), or worldvolume otherwise. We suppose 9Jt to be a D-dimensional vector 
space endowed with the Lorentzian metric ds<xn of signature {1, D — 1); that is 

/ -1 ... \ 

V ••• 1 / 

We consider the toy model of a closed bosonic string. In order to see its dynamics we 
construct dynamical variables x'^((t") where fi = 0,1, . . . , D — 1 and = t, = a 
are coordinates on the canonical cylinder C = {r G (— oo,oo),cr e [0,27r)}. The 
intrinsic metric dsc on C is given by 



We reserve the notation rj^^, for the metric in D-dimensional space and rjap on the 
cylinder C. The dynamical variables describe a C°°-smooth embedding of C in OJl, 
and this allows us to consider the worldsheet (ws) as an embedded manifold in ^M. 
The induced metric ds^g on the worldsheet is given by 

dsls = '^^'•'g^'Q^ da°'da'^ = gap da'^da^. 

The dynamical variables x^(r, a) satisfy the equations of motion derived by making 
use of the least action principle. In Lagrangian mechanics, a system with a configu- 
ration space is characterized by its Lagrangian L, a smooth real-valued functional 
on the direct product of the tangent bundle TOT and the real axis. Let POT stand 
for the path space, that is the space of all paths having the fixed initial point ('ini') 
and the end point ('end'). This space has a structure of an infinite-dimensional 
Frechet manifold. Then, the action functional 5* is a real-valued integral functional 

defined by = Ldt, where t is a real curve parameter. Turning to the spacetime 
configuration space 971 and a worldsheet in 971 we work with the action functional 
S = J^^ ^{x)dx where the Lagrangian L is included in to this spatial integral. So 
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defined action has some pliysical and geometrical relevance. In classical mechanics, 
action, e.g., is the difference between kinetic and potential energies. The minimizing 
curve for the length functional and for the action given by the Lagrangian, which is 
the kinetic energy, coincide. There is no a direct analogue to energy in the relativis- 
tic mechanics. Therefore, Nambu in 1970 |35] and Goto in 1971 {17] proposed to 
choose the area of the worldsheet as the simplest action. This action, known as the 
Nambu-Goto action, admits the form 



Sng = ~T J dx = —T J da"^ \ det gai3 \ 



ws C 

where C is the canonical cylinder and the quantity T, the string tension, has dimension 
mass per unit length. Observe that (-) in front of the integral is due to the Lorentzian 
metric. The string possesses geodesic motion, i.e., the dynamical variables satisfy the 
Euler- Lagrange equation 

1 d f I dx^\ 

where g""^ and gai^ are obtained by rising and lowering of indices with respect to 
the Lorentzian metric. The metric gap is defined on the worldsheet embedded into 
971. It depends on the variables x^, and the equation (12.11) for x^ is highly non- 
linear. Moreover, the square root within the integral results in difficulties in further 
quantization. 

To overcome these problems Polyakov proposed to introduce an analogue of La- 
grange multipliers in 1981 [38] . In fact, he considered a variable metric ds of index 
1 on the worldsheet as on a 2-D manifold, such that ds^ = hapiT, a) da'^da^, and the 
Dirichlet integral 

T f , 2 — -.^.dxf'dx'' 



2 J y^——^^ dcy-daP 

c 

as an action. The Euler-Lagrange equation (regarding to the variation of the dynam- 
ical variables) for the Polyakov action Sp is formally the same as for the Nambu-Goto 
action Sng 

(2.2) , ^ f./|det K3\ /^"^^l = 0, 

but h does not depend on x any longer. So (12. 2p are linear equations for a:;^(r, a). 
Moreover, the Polyakov action is quantized easier due to its linearity. 

There are three degrees of freedom in h because it is symmetric. They can be 
removed by using the equation of motion for h 
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where the left-hand side is a functional derivative. The equations fl2.2H2.3l) perform 
the necessary conditions for minimizing the action. The two-dimensional energy- 
momentum tensor is defined as 

_ -2 1 6Sp 

and the equation of motion implies Tap = 0. Moreover, Sp = Sng under Tajs = 0, 
whereas in general, Sp > Sng- 

A guiding principle in physics is that symmetries in the geometry of spacetime 
define the standard conservation laws. There is no single theorem describing such a 
connection, however this principle is sometimes referred to as the 'Noether theorem', 
although the Noether 1918 theorem [36] itself contains only a very partial statement 
of it. 

By symmetries for the Polyakov action Sp we mean transformations keeping the 
action invariant. They are: 

• Global symmetries 

— Poincare invariance 

x'^ ^ + wf^x'^; 

where cu^ = — cu^ are infinitesimal Lorenzian transformations. 

• Local symmetries 

— reparametrization invariance by diffeomorphisms in 2-D: a"" cr°(r, o") 
implies the invariance of the element 



da^y I det h\ = rfaV|det/i|. 

Weyl rescaling 



Weyl rescaling leaves y^|~det7i[ hap invariant. 

The Poincare and reparametrization invariance of Sp allow us to choose a gauge 
in which the three independent components of h are expressed with just one function, 
typically the conformal flat gauge ha/s = e^'^'^'^^'riap. Substituting this in the Polyakov 
action we obtain 

T f . . dx^ dx^ 



c 



so the gauge fixed action is a quadratic functional with respect to x. Varying it with 
respect to x we arrive at a free wave equation of motion 

(2.4) - (x'^)" = 0, 

where x is the r- and x' is the a- derivative. 
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Weyl invariance implies that Sp does not depend on pij^a). Therefore, the vari- 
ation 6Sp/6p = h^f^Tajs = vanishes that makes Ta/3 traceless. 

Keeping in mind that Ta/3 = one can derive the constraints Tqi = Tio = x^{x^)' = 
and 

(2.5) Too = Tn = ^{x^x^ + x'^^x^') = 0. 

This yields the equations {x ± x')^ = 0, which are called the Virasoro constraints. 
The equation (12. 4p with the constraints (12. 5p . subject to some boundary conditions 
describes the motion of a bosonic string. 

Let us use periodic boundary conditions x^(r, a) = 2;^(r, a + 27r). The general 
solution to the wave equation of motion is 

x^(r, a) = x^(r - a) + <(r + a). 

Let us show that the center of mass of the worldsheet moves as a free particle. Satis- 
fying the boundary condition and separating the linear part, we use the Fourier series 
representation of the solution as 



2 ^ 27iT' V2TTT^n 
x'^(r + a) = -x^ + -^(t + a)p^ + —L= V 



-inir+cr) 



1 



where Xq is the center of mass and is the momentum. The functions and 
are real, and hence, ct^ = a^^ and = [3^^. These coefficients are known under the 
name of oscillators in physics. 

The position of the center of mass of the string is calculated as 

1 f'^'^ 

X^ = - x^[r,a)da = x^, + ^T, 

that shows that the center of mass moves as a free particle starting from Xq. The 
momentum of the center of mass is written 

= / Wda = - / x^da = p^". 







2 

In classical Hamiltonian mechanics the action is a time-integral of the La- 

to 



grangian S = j^^ L{x, x)dt. In our case the Lagrangian becomes 



The Hamiltonian function is given by 

r'27r rji (>2tv 



H = {x>'n^)da -L = - (x^x^ + (x^)'x;)rfa. 
Jo 2 Jo 
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Using an identity for the Kronecker delta 



-/ 



we obtain the Hamiltonian function in terms of oscillators as 



2 

where we set Uq — Pq — -^P^- 

The standard Heisenberg-Poisson-Dirac bracket on phase coordinates looks hke 

{a;^(r, a), a;"(r, a')} = {i:^(r, a), i;"(r, a')} = 0. 

It is convenient to turn to light-cone coordinates on C assuming light speed to 
be 1. They are (± — a ± t and the flat metric ds^ becomes ds^ — d(^dC_. The 

metric components of rjap in the light-cone coordinates are 77++ = 77 = and 

77+- = 77_+ = |. The differential operators become d± — ^{da- ^ 8-^)- 

The Virasoro generators are defined by 



^ -/n ^ ^n, 

27r 



1 /" 1 

= - / r__e^-(-+'^)da = - J] 

nez 

where T_|_+ and T are the principal diagonal components of the energy-momentum 

tensor in the real light-cone coordinates and — (i^ — i^V^ ■ The commutator 
relation for L„ given by the Poisson structure is 

{Lm, Ln\ = i{n — Tn)Ln+m, {Lm, Ln} = — Tn)Ln+m-, {-^mj -^n} = 0. 

The next useful operation is so-called Wick rotation when the non-physical time 
T is replaced by the imaginary time it. Then the light-cone coordinates change to 
C+ — > C_ — > where z — a + it. The differential operators become 9+ — > dz, 
d- dz- Then the Polyakov action in the conformal gauge reads as 



f dz Adz _ „ _ „ 
Sp = -2T I -——rj^^dzx'^dzx'', 



c 

and the equation of motion is Laplacian dl^x^ — 0. The energy- momentum tensor 
has the following components 

Tz, = - ixniix"")' - tx") = Too + 2iT,o, 

T-z-z = Ivf^^iix^y + Unii^n' + ii") = Too - 2irio, 
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-tnz 



and Tzz = T^z = 0. The Virasoro constraints in these coordinates are written as 
Tzz = T-z-z = 0. 

The invariance of Sp with respect to an infinitesimal reparamctrization cr" 
~a _|_ imphes the conservation law for the energy-momentum tensor V^T^is — 
which in ^-coordinates becomes 

d-zTzz + dzT-zz = dzTzz + d-zTzz = 0, 

and making use of Tzz = Tzz — 0, we conclude that Tzz is analytic and Tzz is antian- 
alytic components of T. 

Then in complex coordinates 



V47rr ^ ' 

The Wick rotation leads to the complex coordinates and we conclude that in 
complex coordinates the Virasoro generators L„ are the coefficients of the formal 
Laurent series 

and the commutator relations for L„ in these coordinates are 

{Lrm L„} — i{n — rn)Ln+rm 
so span{L„} forms the Witt algebra. 

2.2. Canonical quantization. In classical mechanics the motion is completely 
defined by the position x and momentum p, which are the elements of a symplectic 
manifold. They specify the state of the classical mechanical system. There is a specific 
function, the Hamiltonian function if, that defines the motion by the equations x = 
VpH, p = —VxH. The symplectic structure (also known as the canonical structure) 
defines the Poisson bracket by {x,p} = 1. 

The main goal of the classical mechanics is to find the trajectories of motion of 
particles, which are the solutions to a Hamiltonian system. The task of quantum me- 
chanics differs because we have no trajectories any longer and a particle is found at a 
point of spacetime with some probability (complex in general). The entire information 
about the quantum system is given by the spectrum and the corresponding collection 
of eigenfunctions for a quantum operator. The initial task of the first (canonical) 
quantization is to find operators X and P corresponding to the dynamical variables 
X and p, which act over an infinite dimensional Hilbert space of smooth functions, 
quantum states. The next step is to find quantum analogue of the Poisson bracket 
and the operator, the Hamiltonian, which defines the total energy of a system, or 
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in other words, an analogue of the Hamiltonian function. Finally, the states of the 
quantum system can be labelled by the eigenvalues of the Hamiltonian. 

Let us start with the definition of the operators X and P. In quantum mechanics 
we replace x ^ -x = X , p ^ —ihVx = P- Acting on quantum states they result in 
Xip = xijj, Pip = —iKVxip- There exist many reasons of this replacement. Let us give 
a toy reason related to the wave function, which appears in the harmonic oscillator 
with the phase = kx — ut, where u is the frequency and k is the coefficient of 
displacement (Hooke's law). The Planck constant is the proportionality constant 
between energy {E) of a photon and the frequency of its associated electromagnetic 
wave: E ^ hv, H ^ 6.62606896(33) x 10'^y27r J s (Joule x seconds). The value Hk 
is the momentum of the displacement in a given direction x. Then the expression 
hdip = hkdx — hjjdt can be interpreted as hdip = pdx — Hdt — dS, where S is the 
action. The real wave function is cos (p, and the complex wave function is z/' = exp(-|S') 
with the amplitude 1. We consider the simplest case in which the amplitude is 
constant=l. Generally, it can be some non-negative function. 

Let us consider the action S, i.e., the time integral of the Lagrangian, on the 
extremal trajectories, the solutions to the Eulcr- Lagrange equation. Then the re- 
sult is a function S{to,t,Xo,x), that depends on the initial and final time to and 
t, and on the initial and final point of the trajectory xq and x. Differentiating 
with respect to the final point x gives p — VxS, where the momentum p is the 
solution of the corresponding Hamiltonian system x = VpH, p = —VxH taken 
at the final time. So we immediately obtain —iKV^ = and the momenta 
represent eigenvalues of the operator —ihV. Thus, the quantum commutator is 
[X, P]ip = {XP - PX)'ip = {x{in)V - ihVx)'ip = ihip, or [X, P] = ih. The quantum 
Poisson brackets are to preserve the classical properties in which {x,p} — 1, therefore 
we come to the conclusion 

quantum 

{^,P}PB- ^ . 

Moreover, in the canonical quantization the fields are replaced by Hermitian op- 
erators. Rescahng our bosonic string, such that ft ^ 1, T ^ 1, we have 

{;r^(r, a), i'^(r, a')} = -ir]^''''6{a - a'), 

{x^iT,a),x''iT,a')} = {,T^(r,a),i;^(r,0} = 0. 

When the canonical quantization procedure is applied in the frames of a given 
quantum field theory, the classical field variable becomes a quantum operator, which 
acts on a quantum state to increase or decrease the number of particles by one. For 
a boson there are two operators: the boson's creation operator and the boson's 
annihilation operator b (commonly known as 'ladder operators') for the system with 
one degree of freedom, or b|. and bk for fields, k E Z. Each operator creates or annihi- 
lates particles in a certain state ip. The next step in quantization is establishing the 
normal ordering operation between creators and annihilators following the Heisenberg 
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uncertainty principle expressed in above Poisson bracket. The commutation relations 
for creators and annihhilators are {b^, b^} = 0, {b, b} = 0, and {b, b^} = 1. The latter 
can be rewritten as bb^ = b^b + 1. A product of creation and annihilation operators 
is said to be in the normal order when all creation operators are to the left of all 
annihilation operators in the product. The process of putting a product into normal 
order is called normal ordering. The normal ordering operation is denoted by : ab :. 
In the case of two boson's operators : b^b := b^b but : bb^ := b^b and bb^ =: bb^ : +1. 

We define the Virasoro generators in the quantum system by introducing normal 
ordering for positively indexed oscillators a„, as annihilators and for the negatively 
indexed ones as creators: 

T -l^T. 



Lm r\ ^ ^ • Pm—nPn 



2 

where m ^ 0. The only operator Lq produces an anomaly because it is already in 
normal order, 



1 1 V- 1 1 V- 

J^O = ^ttottO +92^- «-n«n := ^«OaO + 9 Ck-nttr 



In order to keep the commutation law for oscillators one must subtract from Lq some 
constant c when taking commutator of symmetric operators L„ and L_n- Taking into 
account quantum brackets one arrives at the Virasoro commutation relation 

(2.6) {Lm, Ln}yir = {u - m)Ln+m + ^^'^('^^ " l)5n -m, 

where the constant c is called the central charge and in this particular case it coincides 
with the dimension D. 

The Virasoro constraints can not be written as Ln\ip) = for all n, because 

L.JIV^) = 2n(^|Lo|^) + ^n{n' - ^ 0. 

So the correct Virasoro constraints are imposed by Lnlip) = for n > and (Lq — 
cjltp) = 0. Here {ipl and {ip) are the standard Dirac notations of bra and ket vectors 
of states. Our system is in the physical state \ip) and c = D. 

Early string theory proposed by Yoichiro Nambu and others in 1970 was only 
bosonic. Pierre Ramond, Andre Neveu, and John Schwarz completed the theory by 
inventing fermionic strings to accompany the bosonic ones. 

3. KdV and Virasoro algebra 

The KdV equation appeared in a paper by Diederik Johannes Korteweg and his 
student Gustav de Vries [25j in Philosophical Magazine, 1895, and originally described 
the solitary wave discovered by the Scottish engineer John Scott Russell about half a 
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century earlier. Later on, it became the condition on the potential of the spectral sta- 
bility of the Schrodinger operator, eigenvalues of which represent energy quantization 
for bound systems. The existence of quantized energy levels is verified experimentally 
by observation of the energy emitted or absorbed when the system makes a transition 
from one level to another. 

3.1. Hamiltonian dynamics and integrability. In order to speak on integra- 
bility of KdV let us introduce this notion briefly for finite dimensional Hamiltonian 
systems which will be used in the sequel. There exists a vast amount of modern 
literature dedicated to different approaches and definitions of integrable systems (see, 
e.g., [SSIIZI^)- 

The classical definition of an integrable system in the sense of Liouville applied to 
a Hamiltonian system says, that if we can find independent conserved integrals which 
are pairwise involutory (have vanishing Poisson brackets), this system is partially or 
completely integrable depending on the number of pairwise involutory integrals (see 
e.g., [3], 131, [7] ) . That is each first integral allows us to reduce the order of the system 
not just by one, but by two. We give now the precise definitions. 

Let {N,u!) be a symplectic manifold with the symplectic 2-form u and if be a 

C°°-function defined on N. Then we write H to denote the Hamiltonian vector field 
associated to H. By definition if is a vector field on N, such that 

uj(V,H) = dH(v), 

for every smooth vector field V on A^. In this case the function H is called the 
Hamiltonian function. If H,K E C°°{N), then the Poisson bracket {H,K} is the 

directional derivative of K in the direction of ii, i. e. 

{H, K} = dK(H) = uj(H, K). 

The Jacoby identity and the Leibniz property show that the map H H is a Lie 
algebra homomorphism 

(C-(iV), {•,•}) to (V-(iV), [-,■]), 
where V^lN) denotes the set of vector fields on A^. An equation of the form 

(3.1) x{t) = Hixit)), 

is called the Hamiltonian system. 

Let us consider the particular case of ra-dimensional complex vector space C". 
In this case the cotangent bundle T*(C") for C", that is, isomorphic to C^", has a 
natural symplectic structure. Relatively to the coordinate chart 

T*(n = (zi,...,^„,^i,...,^„), 
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we have the formulas 



n 



'k, 



k=l 




OH OK dH OK 

dz^ dipk dipk dz^ 



(3.2) i'it) = T^, Mt) 



dH 



k = l,2...,n. 



and, using the Poisson brackets, we get 



(3.3) z\t) = {z',H,}, Mt) 



{iJk,H}, k=l,2...,n. 



The latter system has n degrees of freedom. A smooth function $ satisfying {H, $} = 
is called the first integral of the Hamiltonian system. In particular, {H, H} = 0, 
and the Hamiltonian function if is a first integral of the system (13. ip . If the system 
(13.31) has n functionally independent first integrals $ which are pairwise 
involutory, that is $j} = 0, k, j = 1, . . . ,n, then the system is called completely 
integrable in the sense of Liouville. The function H is included in the set of the first 
integrals. The classical theorem of Liouville and Arnold [3J gives a description of 
the motion generated by the completely integrable system (13. 3p . It states that such 
a system admits action-angle coordinates on a connected regular compact invariant 
manifold. 

3.2. From KdV to Virasoro. Let us consider the phase coordinates (field vari- 
ables) defined on the spacetime 5*^ x M, which are from the Sobolev space if" (5^), 
a > —1, regarding to x & for each fixed t G M, i.e., real valued distributions on 
the unit circle u{e^^,t) . Let us simplify u u{x,t), where the new u becomes a 27i 
periodic smooth in x function. Observables are functionals on H°'{S^). 

The famous KdV equation is Ut = 6uu' + u'" on the unit circle can be viewed as 
an integrable Hamiltonian system of infinite dimensions. Several Poisson structures 
can be defined on the space of observables. One was proposed by Gardner [14] . and 
independently, by Zaharov and Faddeev [5T], see others in, e.g., [HI 1181 131j . 

Let us consider two functionals F{u) and G{u), where u G if"(S'^) with the 
vanishing mean value. Expanding u into the Fourier series we obtain 



where U-n = Un, and Z' = Z \ {0}. Then let us consider the functionals F and G as 
functions F = F{qi, g2, • • • P2 • • • ) and G = G{qi, 52, • • • ,Pi,P2 • • • ) with respect 



u 



tnx 
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— Un/n and Pn = U-n-i ^ £ The canonical symplectic form dq A dp implies the 
Poisson structure 



IF G\ = - S^ ( —— = - S^ 



dF dG 

dUn du^n 



Observe that uq can be included now in the latter sum. 

Take into account a useful formula for variational derivatives. If 

f{u, u', u", . . . )dx, 



then 



SF df d df ^ d^ df 



6u du dx du' dx"^ du" 
Hence, if the function u depends on a parameter a, then integrating by parts we 
conclude that 

da Jq 6u da 

In particular, 

OF r^"" SF , 6F 1 ^ OF 1 x-^ dF 



dun Jo 5u^ 5u 27rf-^,dun^ 27r 4^, 



^—mx 



dun Jo Su Su 27r du^ 27r du_n 

Substituting this in the formula for the Poisson bracket we obtain 

IFG\ = — / V n^-^dx = — /" V ^e'nx_d^ -inx^ 



2lT 2-K 

EdF >, dG ^ -ikx^ ■ f ^ d 

rill I ^ rill , I ^ii rlnr ^ii ' 



47^2 / z-^ Qy^ Z-^ du_k J 5u dx Su 



which is true for all functionals of the above integral form. 

If we choose the Hamiltonian function in the form H — —i J^^ {^{u'^) +u^)dx, 
then the evolution equation ii = {u, H} admits the form 

d SH . ... 

u = i—^- — ouu + u . 
dx ou 

One of the most important features of this equation is that it possesses an infinite 
number of conserved quantities (first integrals) e.g., 

/»27r /»27r /'27r -y 

I-i — / udx, Iq = u^dx, h= {-{u'"^) + u^)dx, . . . 

.In .In .In 2 



/.27r 

...,/= / polynomial (^, ■u)dx. 
Jo ^ 
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which are all in involution. One observes the relation of this fact to the Liouville 
integrability for Hamiltonian systems. However, the proof of complete integrability 
is rather difficult task, which was fulfilled in [141 151j . 

If we consider the conserved quantities as Hamiltonians, then we obtain a hierarchy 
constructed as 

u = {u, -.4} = - — , 

which is called the KdV hierarchy. 

The above Poisson structure considered on the phase space formally can be rewrit- 
ten (modulo a constant factor) as 

[u{x),u{y)] =-f6'{x -y), 7 > 0. 

It is degenerate because the center is one-dimensional and the admissible element 

u{x)dx. 



commutes with all observables. Fixing P we get a submanifold which is symplectic. 
The Poisson structure on the phase space gives the Lie structure on the space of 
observables. 

Hamiltonian itself is an observable, and for example, consideration of 

H = — / u^{x, t)dx 
27 Jo 

(corresponding to Jq) gives a trivial evolution equation of motion ii = {H, u} = —u', 
u = u{x — t). The infinite number of evolution equations generated by all integrals is 
the above KdV hierarchy. The Miura transformation s = + u' reduces the modified 
KdV equation ii = u'" + v^u' to the usual KdV and leads to the Hamiltonian 

r.27r 



1 r 

H = — {u^ + u')dx. 
27 Jo 



The Poisson brackets for s satisfy the relation 

{s{x), s{y)} = 7(2(s(x) + s{y)) + 5\x - y) + 5"\x - y). 



which is called Magri brackets 

Fourier coefficients of s{x) are given as 

'.27r 



27r Jo 



The corresponding Lie brackets are {Lm, -^^njvir where c is the central charge, c = ^ 
(J.-L.Gervais [16]). or taking into account quantum effects c = 1 + 6(^ + ^ + 2), see 
[10] (one can observe here the famous gap in CFT for real c G (1, 25)). 
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3.3. Prom Virasoro to KdV. Let us define 

Then, using S{x) = ^Ylin&^^^^ ^^"^ Virasoro commutation relation (12.61) . we 
obtain 

{u{x),u{y)} = —{-6"'{x -y)+ Au(x)6'{x - y) + 2u'6(x - y)). 
c 

Taking /o = | u'^dx, we obtain 

ii = — {u, Jo} = u'" + 6uu'. 

DVT 

KdV as a non-linear PDE is related to the classical and quantum field theories in 
which the infinite number of degrees of freedom follows from the infinite number of 
degrees of freedom for the initial conditions. So it is not surprising to see relations 
between the Virasoro algebra and KdV. As we shall show in forthcoming sections, 
problems of completely different nature (the Lowner-Kufarev evolution) but also of 
the infinite number of degrees of freedom, lead to a rigid algebraic structure given by 
the Virasoro algebra. 

4. Realization on the unit circle 

Mathematically, the Virasoro algebra appeared first as a central extension of the 
Lie algebra of smooth vector fields 0^ on the unit circle (see |15j). Let us denote 
the Lie group of C°° sense preserving diffeomorphisms of the unit circle by Diff S^. 
Each element of Diff 5*^ is represented as 2; = e*"*-^-* with a monotone increasing 
real- valued function a{6), such that a{6 + 2ti) = a{6) + 21^. The space of smooth 
vector fields on a differentiable manifold forms a Lie algebra, where the Lie bracket 
is defined to be the commutator of vector fields. The relation of this Lie algebra to 
Diff is subtile. The Lie algebra to Diff can be associated with the left-invariant 
vector fields Vect S^. But the exponential map, which is the same as the exponential 
map from the tangent space at the origin, is not even locally a homeomorphism. The 
infinitesimal action of Vect 9 ^ 9 + 6(f)(9). To we associate the vector field 
0^, and the Lie brackets are given by 

(4.1) [01, 02] = 0102 - 0201. 

The Virasoro algebra is the unique (up to isomorphism) non-trivial central extension 
of Vect by M given by the Gelfand-Fuchs cocycle |15j . 

4.1. Canonical identification. As an infinite dimensional Lie-Frechet group, 
Diff undergoes certain irregular behaviour, in particular, the exponential map from 
Vect is not a local homeomorphism. The entire necessary background of unitary 
representations of Diff is found in the study of Kirillov's homogeneous Kahlerian 
manifold Diff S^/S^. We deal with the analytic representation of 
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Diff S^/S^. Let S stands for the whole class of univalent functions / in the unit 
disk U normalized by f{z) = z{l + Cn^"") about the origin. By S we denote the 
class of functions from S smooth (C°°) on the boundary 5*^ of U. Given a map / G S 
we construct the adjoint univalent meromorphic map 

g{z) = diz + do + + . . . , 

z 

defined in the exterior U* = {z : \z\ > 1} of U, and such that C\f{U) = g{U*). Both 
functions are extendable onto S^. This conformal welding gives the identification of 
the homogeneous manifold Diff / with the space S: S 3 / o ^[^.i G 

Diff /S^, or with the smooth contours F = f{S^) that enclose univalent domains 
Vt of conformal radius 1 with respect to the origin and such that oo ^ G f2, 
see [I], [23] . Being quasicircles, the smooth contours allow us to embed Diff / 
into the universal Teichmiiller space making use of the above conformal welding, see 
[45] . Coefficients of the univalent functions from S are the natural coordinates on 
the Teichmiiller space. So one can construct complexification of Vect and further 
projection of the holomorphic part to the set M. C which is the projective limit 
of the coefficient bodies = lim„^oo-^n, where 

(4.2) ^„ = {(ci,...,c„): /gS}. 

This construction relates the Kahler structure of both manifolds. The holomorphic 
Virasoro generators can then be realized by the first order differential operators 

oo 

= dj + Y^{k + l)cfca,+fc, J G N, 
fc=i 

in terms of the affine coordinates of M., acting over the set of holomorphic functions, 
where dk = d/dck- We explain the details in the next subsection. 

4.2. Complexification. Let us introduce local coordinates on the manifold 
A4 = Diff S^/S^ in the concordance with the local coordinates on the space S of 
univalent functions smooth on the boundary. Observe that is a real infinite- 
dimensional manifold, whereas S is a complex manifold. We are aimed at a com- 
plexification of TAi which admits a holomorphic projection to TS, where as usual, 
Vect o'S'^ = Vect S^/const is a module over the ring of smooth functions, which is 
associated with the tangent bundle TA4. 

Two operations are to be considered: complexification, conjugation, and almost 
complex structure. Given a real vector space V the complexification Vc is defined 
as the tensor product with the complex numbers V (S>m C, that often is written as 
Vc = VQ)iV. The subscript M indicates that we take the real tensor product, we omit 
it in the sequel. Elements of Vc are of the form v^z. In addition, the vector space Vc 
is a complex vector space that follows by defining multiplication by complex numbers, 
a{v^z) = v®az for complex a and z and v & V. The space V is naturally embedded 
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into V" (8> C by identifying V with V ^ 1. Conjugation is defined by introducing a 
canonical conjugation map on Vc as v ^ z = v ® z. 

An almost complex structure J on K is a linear transformation J : V ^ V such 
that J2 = -/. It can be extended by linearity to the complex structure J on Vc by 
J{v ® z) = J{v) ® z. Observe that 



J{v ® z) = Jv ® z = Jv ® z = J{v ® z) = J{v ® z). 

Eigenvalues of extended J are ±z, and there are two eigenspaces and 
corresponding to them given by projecting ^(1 =F iJ)v. Vc is decomposed into the 
direct sum Vc = © \/(o.i)^ where = {v ® 1 - J{y) ® i\v G V^} and 

y(o,i) _ 1^ (g, ]^ _|_ j^^^ (g)z|t> G V}. In the case of existence of such a complex structure 
J, the vector spaces V^"^'^^ and V^'^'^^ give complex coordinates on V . 

An almost complex structure on Vect o'S'^ may be defined as follows (see [1]). We 
identify Vect o'S'^ with the functions with vanishing mean value over S^. It gives 

oo 

(j){9) = ttn COS nO + bn sin n9. 

n=l 

Let us define an almost complex structure by the operator 

oo 

J{(j)){d) = — a^isin n6 + 6„ cos n6. 

n=l 

Then = —id. On Vect o^*^ ® C, the operator J diagonalizes and we have the 
identification 

^ oo 

Vect oS^3(P^v:= -(0 - iJ(0)) = ^(a„ - i6„)e*"^ G (Vect ® C)(^'°\ 

n=l 

and the latter extends into the unit disk as a holomorphic function. 

The Kirillov infinitesimal action [22] of Vect on S is given by a variational 
formula due to Schaeffer and Spencer [44\ page 32] which lifts the actions from the 
Lie algebra Vect onto S. Let / G S and let 0(e*^) := (j){e) G Vect be a 
real- valued function in ^ G (0,27r]. The infinitesimal action 9 9 + ecpi^e^^) yields a 
variation of the univalent function f*{z) = / + e5yf{z) + o(e), where 

(4.3) 5J{z)^f"^'^ f f^f'HV -Hd^ 



27r i V /H ; w{f{w)-f{z)y 

and f by the above identification. Kirillov and Yuriev 23]. |24j (see also [1]) 
established that the variations S^f{() are closed with respect to the commutator (14. ip . 
and the induced Lie algebra is the same as Vect O'S'^. The Schaeffer- Spencer operator 
is linear. 

Treating TAi as a real vector space, the operator 5,^ transfers the complex struc- 
ture J from Vect o'S'^ to TM. by J{S^) := ^j(</,). By abuse of notation, we denote 
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the new complex structure on TAi by the same character J. Then it sphts the 
complexification TA4c into two eigenspaces TA4c = TAi^^'^^ (B TAi^'^'^\ Therefore, 
= 5(f,_ij{^^) := 5^ — iJ{5^) G TJ^^^'^\ Observe that 2zdz = —idg on the unit 
circle z = e*^, and Lk = z^^^d/dz = —^ie'^''^d/d9 on S^. Let us take the basis of 
Vect qS^^C in the form = —ie^^^ in order to keep the index of vector fields the same 
as for Lfc. Then, the commutator satisfies the Witt relation {z/^, Vn} = {n — rnjVn+m- 
Taking elements z/^ = —iw^, |w| = 1 in the integrand of (14.31) we calculate the residue 
in (14.31) and obtain so called Kirillov operators 

L^[f]{z) = 6,J{z) = z^^'f{z), J = 1,2,..., 

which are the holomorphic coordinates on TA4^^'^\ In terms of the affine coordinates 
in Ai we get the Kirillov operators as 

oo 

= + + l)cfc5j+fc, 

k=l 

where = d/dck. They satisfy the Witt commutation relation 

For k = we obtain the operator Lq, which corresponds to the constant vectors 
from Vect S^, Lo[/](-2) = zf'{z) — f{z). The elements of the Fourier basis —ie~'''^^ 
with negative indices (corresponding to TAd'^'^'^^) are extended into U by —iz^^. 
Substituting them in (14.31) we get very complex formulas for which functionally 
depend on (see [T], |22] ) and might play the role of conjugates to L^. The first 
two operators are calculated as 

L_^[f]{z) = /'(z) - 2ci/(^) - 1, 

L-2[f]{z) = ^^--^-3ci + (c?-4c2)/(z), 

see [24] . 

This procedure gives a nice links between representations of the Virasoro alge- 
bra and the theory of univalent functions. The Lowner-Kufarev equations proved 
to be a powerful tool to work with univalent functions (the famous Biberbach con- 
jecture was proved [5] using Lowner method). In the following section we show how 
Lowner-Kufarev equations can be used in a representation of the Virasoro algebra. In 
particular, we identify TAi^^''^'' with TAi, equipped with its natural complex structure 
given by coefficients of univalent functions, by means the Lowner-Kufarev PDE. 

5. Lowner-Kufarev Equations 

A time-parameter family Q{t) of simply connected hyperbolic univalent domains 
forms a Lowner subordination chain in the complex plane C, for < t < r (where r 
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may be oo), if Q{t) ^ ^{s), whenever t < s. We suppose that the origin is an interior 
point of the Caratheodory kernel of {^^(t)}[=o- 

A Lowner subordination chain Q{t) is described by a time-dependent family of 
conformal maps z = f{(, t) from the unit disk U = {C, : |C| < 1} onto normalized 
by /(C,t) = ai(t)C + a2{t)C + • • • , ai{t) > 0, ai(t) > 0. After Lowner's 1923 seminal 
paper |27j a fundamental contribution to the theory of Lowner chains was made by 
Pommerenke [391 140] who described governing evolution equations in partial and 
ordinary derivatives, known now as the Lowner- Kuf are v equations due to Kufarev's 
work [26j. 

One can normalize the growth of evolution of a subordination chain by the con- 
formal radius of Q{t) with respect to the origin by ai(t) = e*. 

Lowner [27] studied a time-parameter semigroup of conformal one-slit maps of 
the unit disk U arriving then at an evolution equation called after him. His main 
achievement was an infinitesimal description of the semi-flow of such maps by the 
Schwarz kernel that led him to the Lowner equation. This crucial result was then 
generalized in several ways (see [40] and the references therein). 

We say that the function p is from the Caratheodory class if it is analytic in U, 
normalized as p{C) = 1 -f piC + P2C^ -l- . . . , ( E U, and such that Re p{() > 
in U. Pommerenke [39|. 140] proved that given a subordination chain of domains Q{t) 
defined for t G [0,r), there exists a function p{C,t), measurable in t G [0,r) for any 
fixed z E U, and from the Caratheodory class for almost all t G [0, r), such that the 
conformal mapping f : U ^ Q{t) solves the equation 

(5.1) ^.c^?f^„c,o, 

for C G f/ and for almost all t G [0,r). The equation (15.11) is called the Lowner- 
Kufarev equation due to two seminal papers: by Lowner [27] who considered the 
case when 

(5-2) PiCt) = -^^, 

where u{t) is a continuous function regarding to t G [0,r), and by Kufarev [26j who 
proved differentiability of / in t for all ( from the kernel of {Q{t)} in the case of 
general p in the Caratheodory class. 

Let us consider a reverse process. We are given an initial domain Q{0) = Qq 
(and therefore, the initial mapping /(C^O) = /o(C))5 ^^^^ ^ function p{C,t) of positive 
real part normalized by p{(, t) = 1 + pi( + . . . . Let us solve the equation (15. ip and 
ask ourselves, whether the solution f{C:t) defines a subordination chain of simply 
connected univalent domains f{U,t). The initial condition /(C,0) = /o(C) is not 
given on the characteristics of the partial differential equation (15. ip . hence the solution 
exists and is unique but not necessarily univalent. Assuming s as a parameter along 
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the characteristics we have 

dt dC , . df 

as as as 

with the initial conditions t(0) = 0, ({0) = z, f{(,0) = fo{C), where z is in U. 
Obviously, t = s. Observe that the domain of ( is the entire unit disk. However, the 
solutions to the second equation of the characteristic system range within the unit 
disk but do not fill it. Therefore, introducing another letter w (in order to distinguish 
the function w{z,t) from the variable () we arrive at the Cauchy problem for the 
Lowner-Kufarev equation in ordinary derivatives 

dw 

(5.3) — = -wp{w,t), 

for a function ( = w{z, t) with the initial condition w{z, 0) = z. The equation (15.31) is a 
non-trivial characteristic equation for (15. ip . Unfortunately, this approach requires the 
extension of fo{w~^{(,t)) into the whole U {w~^ means the inverse function) because 
the solution to (15.11) is the function f{(, t) given as fo{ui~^{(, t)), where ( = w{z, s) is 
a solution of the initial value problem for the characteristic equation (15. 3p that maps 
U into U. Therefore, the solution of the initial value problem for the equation (15.10 
may be non- univalent. 

Solutions to the equation (15.30 are holomorphic univalent functions w{z,t) = 
e~^z + 02(^)2;^ + .. . in the unit disk that map U into itself. Every function / from 
the class S can be represented by the limit 

(5.4) f{z) = lim e^w{z,t), 

where w{z,t) is a solution to (15. 3p with some function p{z,t) of positive real part for 
almost alH > (see |40l pages 159-163]). Each function p{z,t) generates a unique 
function from the class S. The reciprocal statement is not true. In general, a function 
/ e S can be obtained using different functions p{-,t). 

Now we are ready to formulate the condition of univalence of the solution to the 
equation (15. ip . which can be obtained by combination of known results of [40j . 

Theorem 1. [401 141] Given a function p{C,t) of positive real part normalized 
by p{C,t) = 1 + piC + • • • ; the solution to the equation ( (5. 1\) is unique, analytic and 
univalent with respect to ( for almost all t > 0, if and only if, the initial condition 
/o(C) is taken in the form (15. 4p . where the function w{(,t) is the solution to the 
equation (15.31) with the same driving function p. 

Recently, we started to look at Lowner-Kufarev equations from the point of view 
of motion in the space of univalent functions where Hamiltonian and Lagrangian for- 
malisms play a central role (see, [46]). Some connections with the Virasoro algebra 
were also observed in |33L 146] . The present paper generalizes these attempts and 
gives their closed form. The main conclusion is that the Lowner-Kufarev equations 
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are naturally linked to the holomorphic part of the Virasoro algebra. Taking holo- 
morphic Virasoro generators L„ as a basis of the tangent space to the coefficient 
body for univalent functions at a fixed point, we see that the driving function in the 
Lowner-Kufarev theory generates generalized moments for motions within the space of 
univalent functions. Its norm represents the energy of this motion. The holomorphic 
Virasoro generators in their co-tangent form will become conserved quantities of the 
Lowner-Kufarev ODE. The Lowner-Kufarev PDE becomes a transition formula from 
the affine basis to Kirillov's basis of the holomorphic part of the complexified tangent 
space to at any point. Finally, we propose to study an alternate Lowner-Kufarev 
evolution instead of subordination. 



6. Witt algebra and the classical Lovi^ner-Kufarev equations 

In the following subsections we reveal the structural role of the Witt algebra as 
a background of the classical Lowner-Kufarev contour evolution. As we see further, 
the conformal anomaly and the Virasoro algebra appear as a quantum or stochastic 
effect in SLE. 

6.1. Lowner-Kufarev ODE. Let us consider the functions 

/ oo 

w{z, t) = e~^z 1 1 + c„(t);z" 

n=l 



satisfying the Lowner-Kufarev ODE 

dw 

(6.1) — = -wp{w,t), 

with the initial condition w{z, 0) = z, and with the function p{z, t) = 1 +pi(t)z + . . . 
which is holomorphic in U and measurable with respect to t G [0, oo), such that 
Re p > in [/. The function w{z,t) is univalent and maps U into U. 

Lemma 1. Let the function w{z,t) be a solution to the Cauchy problem for the 
equation (16.11) with the initial condition w{z,0) = z. If the driving function p{-,t), 
being from the Caratheodory class for almost all t > 0, is smooth in the closure U 
of the unit disk U and summable with respect to t, then the boundaries of the domains 
B{t) = w{U, t) G U are smooth for all t. 

Proof. Observe that the continuous and differentiable dependence of the solution 
of a differential equation x = F{t,x) on the initial condition x(0) = xq is a classical 
problem. One can refer, e.g., to [49] in order to assure that summability of F{-,x) 
regarding to t for each fixed x and continuous differentiability {C^ with respect to 
X for almost all t) imply that the solution x(t,xo) exists, is unique, and is with 
respect to Xq. In our case, the solution to (16.11) exists, is unique and analytic in U, 
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and, moreover, on its boundary S^. Let us differentiate (16.11) inside the unit disk 
U with respect to z and write 

t 

logw' = - {p{w{z, r), r) + w{z, t)p'{w{z, r), r))cir. 



choosing the branch of the logarithm such as logw'{0,t) = —t. This equahty is 
extendable onto 5*^ because the right-hand side is, and therefore, w' is and w is 
on S^. We continue analogously and write the formula 

t 

w" = —w' J {2w'{z, t)p\w{z, t),t) + w{z, r)w'{z, r)p"(w{z, r), T))dT, 



which guarantees that w is on S^. Finally, we come to the conclusion that w is 
on S\ □ 

Let f{z,t) denote e*w{z,t). The limit limt^oo f{z,t) is known [40] to be a repre- 
sentation of all univalent functions. 

Let the driving term p{z, t) in the Lowner-Kufarev ODE be from the Caratheodory 
class for almost all t > 0, C°° smooth in If, and summable with respect to t. Then 
the domains Q{t) = w{U,t) have smooth boundary dfl{t). So the Lowner equation 
can be extended onto the closed unit disk U = U U S^. 

Consider the Hamiltonian given by 

(6.2) H= [ f{z,t)il-pie-'fiz,t),tmz,t)-, 

J 

on the unit circle z & S^, where %l){z,t) is a formal series 

oo 
n=—k 

defined about the unit circle for any k > 0. The Poisson structure on the space 
(/, ip) is given by the canonical brackets 

or in coordinate form (only ipn for n > 1 are independent co-vectors corresponding to 
the tangent vectors dn with respect to the canonical Hermitean product for analytic 
functions) 

dp dq dp dq 



{p, (i} = Yl 



^ dCn dlpn dlpn ' 
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Here 

P(t)= / p{z,t)-, Q{t)= [ qiz,t)-. 
J tz J tz 

The Hamiltonian system becomes 

(6.3) *M = /(I -p(e-/.«)) = ^| = {/.//}. 

for the position coordinates and 

(6.4) ^ = -(1 - t) - e-*/p'(e-V, m = ^ = bP. H}, 

for the momenta, where ^ and 4t are the variational derivatives. So the phase 

coordinates (/, tp) play the role of the canonical Hamiltonian pair. 

The coefficients Cn are the complex local coordinates on A^, so in these coordinates 
we have 

f dz 

Cn = Cn- — I w{z,t)p{w{z,t),t)- 



- f ye-''\e'w)'+%^, n>l. 

•J 1 1 



2m 

s 



1 k=l 



Let us fix some n and project the infinite dimensional Hamiltonian system on an n- 
dimensional Ain- Momenta in coordinates form an adjoint vector 

ipit) = (V^i(t), . . . , ipnit))'^, with complex-valued coordinates ipi, ■ ■ ■ yi'n for any fixed 
n. The dynamical equations for momenta governed by the Hamiltonian function (16.21) 
are 

k=i 

and 

(6.5) $n = 0. 



In particular, 



ci = -e Vi, 



C2 = -2e ViCi - e ^^2, 

C3 = -e"Vi(2c2 + cl) - 3e~^V2Ci - e'^Vs, 
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for n = 3 we have 

ipi = 2e"*piV^2 + (2e"*]9iCi + 3e~^V2)^3, 
4 = 0. 

Let us set the functfon L{z) := f {z , t)ip {z , t) . Let (L(2;))<o mean the part of the 
Laurent series for L{z) with negative powers of z, 

1 - - 1 °°L 

(L(2;))<o = (^1 + 2ci^2 + Sc2iJ3 + ...)- + {^2 + 2c,tlj3 + •••)- + ■■• = 5Z 4- 

k=l 

Then, the functions L{z) and (L(z))<o are time-independent for all z G S*^. 

It is easily seen that, passing from the cotangent vectors ipk to the tangent vectors 
the coefficients of {L{z))^q defined on the tangent bundle TA4.^^'^^ are exactly 
the Kirillov vector fields Lk. Corresponding co- vector fields Lk are conserved by the 
Lowner-Kufarev ODE because = {Lk, H} = 0. The above Poisson structure 
coincides with that given by the Witt brackets introduced for Lk previously. For 
finite-dimensional grades this result was obtained in [33] . 
Let us formulate the result as a theorem. 

Theorem 2. Let the driving term p{z,t) in the Lowner-Kufarev ODE be from 
the Caratheodory class for almost all t > 0, C°° smooth in U, and summahle with 
respect to t. Then the Kirillov co-vector fields are the conserved quantities for the 
Hamiltonian system 1^6. 3\ - 6.4 ) generated by the Lowner-Kufarev ODE. 



Remark 1. Another way to construct a Hamiltonian system could be based on 
the symplectic structure given by the Kdhlerian form on Diff S-^/S-^. However, there 
is no explicit expression for such form in terms of functions / G S. Moreover, there 
must be a Hamiltonian formulation in which the Lowner-Kufarev equation becomes 
an evolution equation. This remains an open problem. 

Remark 2. At a first glance the situation with an ODE with a parameter is quite 
simple. Indeed, if we solve an equation of type f{t,e^^) = F{f{t,e^^),t), then fixing 
9 we have an integral of motion C = I{f{t,-),t) = const. Then, releasing 9, we 
have C(e*^) = I{f{t,e'^^),t). Expanding C(e*^) into the Fourier series, we obtain an 
infinite number of conserved quantities, but they do not manifest an infinite number of 
degrees of freedom that govern the motion as in the field theory where the governing 
equations are PDE. In our case, we have not only one trajectory fixing the initial 
condition but a pensil of trajectories because our equation has an infinite number 
of control parameters, the Taylor coefficients of the function p{z,t), which form a 
bounded non-linear set of admissible controls. Therefore, we operate with sections of 
the tangent and co-tangent bundles to the inifinite dimensional manifold M. instead 
of vector fields along one trajectory as in usual ODE. 
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Remark 3. No linear combinations L\ of Li, . . . , Ln, . . . allows us to reduce the 
system of {L^} to a new system of involutory {L*^^} in order to claim the Liouville 
integrahility of our system. Observe that the coefficients in these linear combinations 
must be constants to keep conservation laws. 

6.2. Construction of Lq and Consider again tlie generating function 

L{z) = f'{z,t)ip{z,t) and the 'non-negative' part {L{z))>q of the Laurent series for 

L{z), 

{L{z))>o = + 2ciV^i + ?>C2i)2 + ...) + (^-1 + 2ci^/^o + ScsV^i + ...)z + ... 

oo 



k=0 



All L-k are conserved by the construction. Define ipQ — — Yl'^=i '^k'^k, and 

Lq^ Co - (V'o - V'o)- 

The operator Lq acts on the class S by LQ[f]{z) = zf'{z) — f{z). Next define L_i = 
£_! - (V^-i - Ci) - 2ci(^o - V'o). where ^ = 0. Then, 

L_,[f]{z)^f'{z)-2c,f{z)-l 

Finally, 

L-2 = - (^-2 - ^-2) - 2ci(^_i - - 3c2(V^o - V'o)- 
We choose 'il)*_2 = (ca — 3ciC2 + cf)ipi + . . . , so that 

1 1 



and 



r-2[m = ; - Tf^y - ^1 - - 

L-2[f]{z) = ^ - 7^ - + (c? - 4c2)/(^). 



An important fact is that 

Lq = Cii/Si + 2C2V'2 + ■ ■ ■ , 

L_i = (3c2 - 2c?)^i + . . . , 
L_2 = (5c3 - 6cic2 + 2cl)^i + 
are linear with respect to ■0^, k > 1, and therefore, are co- vectors. Equivalently 

.k _ df 



Lo_i_2[/](^) =function(ci,C2,...)^ + ..., z' 



dck-i 

All other co-vectors we construct by our Poisson brackets as 

1 1 

L-n — — —^{L-n+l, L-i} — — — — {L_„+2, -^-2}- 

The form of the Poisson brackets guarantees us that all L_„ are linear with respect 
to ■01, ■02, ■ ■ ■ and span the anti-holomorphic part of the co-tangent bundle T^^'^) M.. 
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Let us summarize the above in the following conclusion. We considered a non- 
linear contour dynamics given by the Lowner-Kufarev equation. It turned out to be 
underlined by an algebraic structure, namely, by the Witt algebra spanned by the 
Virasoro generators L„, n G Z. 

• Ln, n = 1,2,... are the holomorphic Virasoro generators. They span a 
holomorphic vector bundle over the space of univalent functions S, smooth 
on the boundary. In their co-vector form, L„ are conserved by the Lowner- 
Kufarev evolution. 

• Lq is the central element. 

• L-n, n = 1,2, . . . are the antiholomorphic Virasoro generators. They span 
the antiholomorphic vector bundle. In their co-vector form, L_n contain a 
conserved part and we give an iterative method to obtain all of them based 
on the Poisson structure of the Lowner-Kufarev evolution. 

6.3. Lowner-Kufarev PDE. The Lowner equation in partial derivatives is 
w{C, t) = Cw'iC, t)p{C, t). Re p(C, t) > 0, Id < 1. 

with some initial condition w{z, 0) = fo{z)- Let us consider the one-parameter family 
of functions f{z,t) = e~^w{z,t) = z{l + Yl'^=i Cnit) z"') , f{z,0) = fo{z) as a path 
in S. At the identity id we have that Tj^S = Tj^TVl^^'^^ = Tij^M.. A path in the 
coefficient body M. in the neighbourhood of the identity is (ci(t), . . . , c„(t), . . . ) with 
the velocity vector Cidi -f ■ ■ ■ -|- c„9„ + ■ ■ ■ G Ti^Ai. 

Taking the Virasoro generators {Lk}, k > 1, as an algebraic basis in TA^j-^'°^ we 
wish the velocity vector written in this new basis to be 

(6.6) Cidi H h Cndn H = UiLi + . . . UnLn + . . .. 

We compare (16.61) with the Lowner-Kufarev equation 

(6.7) / = Cidi H h Cndn H = Zf'p{z, t) - f = Lq + UiLi + ... UnLn + . . ., 

where p{z, t) = 1 + Uiz -!-■■■ + m„2;" + . . . , and Lq/ = zf — /. In view of similarity 
between these two expressions (16. 6p and (16.71) . we notice that 

• a new term Lq appears in the Lowner-Kufarev equation; 

• the function p{z, t) with positive real part corresponds to subordination, 
whereas for generic trajectories it may have real part of arbitrary sign. We 
call this an alternate Lowner-Kufarev evolution; 

• the vector Lq corresponds exactly to the rotation: 

e-/(e-z) = f{z) - te{zf\z) - f{z)) + o{e). 

Let us consider the set Sq of non-normalized smooth univalent functions of the 
form F{z, t) = ao(t)z + ai(t)z'^ + . . . , with a tangent vector ao<9o + ■ ■ ■ -|- dndn + . . . , 
where dk = d/dak, k = 0,1,2,.... Our aim is to define two distributions of co- 
dimension 1 for the tangent bundle TSq, that form the tangent bundle TS. This will 
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be realized by means of formulas dSSD and Notice that dkF = z''+\ Setting 

Lfc(F) := z'^+'^F' we get 

F = dodo H h a„<9„ H = zf'p{z, t) = UqLo + UiLi + ... m„L„ + . . . , 

where t) = uq + uiz + ■ ■ ■ + UnZ^ + . . . . This alternate Lowner-Kufarev equation 
represents recalculation of the tangent vector in the new basis 

aodo H h cindn H = UoLq + . . . m„L„ + . . ., 

where = aodk + 2aidk+i + . . . . 

Let us present the distributions. We start with F G So, then we define / G S. 
The necessary distribution is the map 

So3F^TfS^ TpSo. 

The analytic form of the first distribution is the following factorization fi{z,t) = 
—F(z, t) = z + ^z^ + . . . , so that 

(6.8) f^ = ,f[p^z,t)-^f„ 

ao 

where uq = —. Then we obtain 



Cidi H h Cndn H = -^^0 + UiLi H h UnLn + • • • 

where Lq/i = Uoizf[ - fi), Lfc/i = z''+'^f[, Ck = ^, dk = In particular, oq = e* 
implies the Lowner-Kufarev equation for arbitrary sign of Re p. 

The analytic form of the second distribution becomes f2{z,t) = F{—z,t) = z + 



' ao 

-fz-^..., - 



^z'^ + . . . , so that 



(6.9) f^ = zf^p{±,t)-^zf^, 

ao ao 

where again uo = In the coefficient form we get 

Cidi + h Cndn H = UiLi H h Un.Ln + . ■ ■ 

where L^/a = z^^'f^, cj, = dk = ^. 

Observe that the equation (16.91) gives another identification of TAi^^'^^ with TA4. 
Finally, let us make an explicit calculation of Lq, which for ao = e* we continue 
to denote by Lq. Using Kirilov's basis Li, L2, ... as a linear combination we write 



m=l 
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The coefficients 11^ are polynomials, which can be obtained using the following re- 
current formulas 

m— 1 m 
i^l = 0, Km = - ^ Jim -j + l)Cm-jCj, Um = mCm + ^ Km^j+iPj_i, 

where Pk are polynomials 

k 

(6.10) Po = l, A = -2ci, P2 = 4cl-3C2, = - J](j + l)c,P,_,, 

i=i 

Let us summarize the above considerations in the following theorem. 

Theorem 3. The L owner- Kufarev PDF (16. 8p gives the distribution of co- dimension 
1 inside the tangent bundle TSq of non-normalized smooth univalent functions So, that 
forms the tangent bundle TS. 

The equation (16. 9p gives another distribution, and moreover, it makes the explicit 
correspondence between the natural complex structure o/TS and the complex structure 
ofTM^^'^^ at each point / G S. 

One of the reason to consider the alternate Lowner-Kufarev PDE is the regularized 
canonical Brownian motion on smooth Jordan curves. For all Sobolev metrics H2+^^ 
the classical theory of stochastic flows allows to construct Brownian motions on 
diffeomorphism group of S^. The case 3/ 2 is critical. Malliavin [32] constructed the 
canonical Brownian motion on the Lie algebra Vect for the Sobolev norm 

^3/2 _ 

Another construction was proposed in |12j . Airault and Ren [2] proved that the 
infinitesimal version of the Brownian flow is Holder continuous with any exponent 
P<1. 

The regularized canonical Brownian motion on Diff is a stochastic flow on 
associated to the Ito stochastic differential equation 

Cxti^) = , , {x2„{t) COS n9 - X2„-i{t) sin n9), 
, yn^ — n 

n=l * 

where {x^} is a sequence of independent real-valued Brownian motions and r G 
(0, 1) and the series for Cx t(^) ^ Gaussian trigonometric series. Kunita's theory of 
stochastic flows asserts that the mapping 6 ^{6) is a C°° diffeomorphism and the 
limit lim g!^ j = g^^t exists uniformly in 6. The random homeomorphism g^^t is called 

canonical Brownian motion on Diff S^, see [2], I12L 1321 142] . It was shown in [2|, 112] . 
that this random homeomorphism is Holder continuous. 

The canonical Brownian motion can be defined not only on Diff S^, but also on 
the space of C°°-smooth Jordan curves by conformal welding. This leads to dynamics 
of random loops which are not subordinated. 



CONFORMAL FIELD THEORY AND LOWNER-KUFAREV EVOLUTION 



29 



7. Elliptic operators over the coefficient body 

The Kirillov first order differential operators generate the elliptic operator 
^ |Lfcp. In this section we construct the geodesic equation and find geodesies with 
constant velocity coordinates in the field of this operator. In particular, we shall 
prove that the norm of the driving function in the Lowner-Kufarev theory gives the 
minimal energy of the motion in this field. 

7.1. Dynamics within the coefficient body. This subsection is auxiliary. 
Let us recall the geometry of the coefficient body Ain for finite n with respect to 
Kirillov's basis L^, = 1, . . . ,n of the tangent bundle TAin- The affine coordinates 
are introduced by projecting 

oo 

M3 f = z(^l + '^Ckz''^ ^ (Ci,. . .,C„) G Mn. 
k=l 

The manifold A1„ was studied actively in the middle of the last century, see e.g., 
[5], [44] . We compile some important properties of A4n below: 

(i) A4n is homeomorphic to a (2n — 2)-dimensional ball and its boundary dAin 
is homeomorphic to a (2n — 3)-dimensional sphere; 

(ii) every point x G dAin corresponds to exactly one function / G S which is 
called a boundary function for A^„; 

(iii) boundary functions map the unit disk U onto the complex plane C minus 
piecewise analytic Jordan arcs forming a tree with a root at infinity and 
having at most n tips, 

(iv) with the exception for a set of smaller dimension, at every point x G dAin 
there exists a normal vector satisfying the Lipschitz condition; 

(v) there exists a connected open set Xi on dA4n, such that the boundary dAin 
is an analytic hypersurface at every point of Xi. The points of dAin cor- 
responding to the functions that give the extremum to a linear functional 
belong to the closure of Xi. 

Properties (ii) and (iii) imply that the functions from S deliver interior points of 
Ain- The Kirillov operators Lj restricted onto A4n give truncated vector fields 

n-j 

= + + l)Cfc<9j+fc, 

fc=l 

which we, if it causes no confusion, continue denoting by Lj in this section. In [33j 
based on the Lowner-Kufarev representation, we showed that these Lj can be obtained 
from a partially integrable Hamiltonian system for the coefficients in which the first 
integrals coincide with Lj. 

Let c{t) = (ci(t), . . . , c„(t)) be a smooth trajectory in A^„; that is a map 
c : [0, 1] Ain- Then the velocity vector c(t) written in the affine basis as c(t) = 
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Ci(t) di + . . . + Cn{t) d n can be also represented in the basis of vector fields Li^ • • • i Ln 
(compare with (I6.9p ) as 

(7.1) c{t) =ci(t)9i + ... + c„(t)9„ 

= UiLi + + . . . + UnLn, 

where the coefficients Uk can be written in the recurrent form as 

fc-i 

(7.2) Ml = ci, Uk = Ck- ^(j + l)cjUk-j. 

i=i 

Expressing Uk in terms of Ck and Ck, we get 

k-l 

(7.3) Uk = Ck + ^ PjCk-j- 

i=i 

One may notice that these polynomials are the first coefficients of the holomorphic 
function l/f'{z), where / G S. In the infinite dimensional case this follows from 
the Lowner-Kufarev equation (16. 9p with = e*. Kirillov's fields act over these 
polynomials as 

LkPn = {n — 2k — l)Pn^k n > k and L^Pn = n < k. 

Proposition 1. Let us give the conjugate to {Li,...,Ln} basis of one- forms. 
We define 

uJi = dci, 

UJ2 = dc2 — 2ciUJi, 

5 

(7.4) UJn = dCn - ^(j + l)CjUJn^j. 

i=i 

Then 

UJn{Ln) = 1, i^niLk) =0 if k ^ 71. 

Proof. If > n, then the vector fields Lk do not contain dn for k > n. Since 
the form uin depends only on dcj with j < n, then 

n-l 

Un{Lk) = dn{Lk) - ^(j + l)cjUn-j{Lk) = for k> n> n- j. 
i=i 

If n = fc, then 

n-l 

UJn{Ln) = dn{Ln) - ^(j + l)Cjl^„_j (L„) = 1 + for U > U - j. 
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To prove the case k < n we apply the induction. Let us show for Li. We have 

U2iLi) = dc2{Li) - 2ci(Li) = 2ci - 2ci = 0. 
We suppose that LJn{Li) = 0. Then 

n 

a;„+i(Li) = dcn+iiLi) - ^(j + l)cjUJn+i-jiLi) = (n + l)c„ - (n + l)c„u;i(Li) = 0. 

The same arguments work for ujn{Lk) with k < n. □ 

In the affine basis the forms can be written making use of the polynomials Pn- 
We observe that one-forms Uk are defined in a similar way as the coordinates Uk with 
respect to the Kirillov vector fields L^. Thus, if we develop the recurrent relations (17.40 
and collect the terms with dcn we get 

fc-i 

i^fc = dck + ^ Pjdck-j. k = 1, . . . , n. 
i=i 

By the duality of tangent and co-tangent bundles the information about the mo- 
tion is encoded by these one-forms. 

7.2. Hamiltonian equations. There exists an Hermitian form on TAin, such 
that the system {Li, . . . , L„} is orthonormal with respect to this form. The operator 
L = l-^fcP is elliptic, and we write the Hamilton function H{c,c,ip,ip) defined on 
the co-tangent bundle, corresponding to the operator L as H{c, c,ip,'4') = Ylk=i I^^P' 
where 

n—k 

h = i'k + ^(j + i)cj4}k+j. 

The corresponding Hamiltonian system admits the form 

dH - 

Ci = = li 



dH 



L + + ^)Cjln-j 



dH 

A = -j^ = -{p + ^)^hipk+p 

k=i 



t OH 
oc„. 
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Let us observe that 



(7.5) h = $^(j - k)ljl^ 



n—k 

i+k- 



Expressing Ik from the first n Hamilton equations we get 

fe-i 

(7.6) Ik = Ck + ^ PjCk-j, k = l,...,n. 

i=i 

We can decouple the Hamiltonian system making use of (17. 5p and (17. 6p which leads 
us to the following non-linear differential equations of the second order 

k-l k-1 

Ck = lk + ^(j + l)Cjlk-l + ^(j + l)c/fc-Z, 
j=l J=l 

where Ik are expressed in terms of the product of Ijlj+k by (17. Sp . and the last products 
depend on Pj, Pj and c, Cj for the corresponding indices j by (17. 6p . For example, 



n-l j-1 



Cl 



ii = ^(j - 1) (^cj + J2 Pp^j~p) (cj+i + 

j = l p=l q=l 



1 / ■ 



Comparing (17. 6p and (17. 3p . we conclude that Ik = Uk and Uk satisfy the differential 
equations 



(7.7) ilk = - k)ujUj 



n—k 

+k, 



on the solution of the Hamiltonian system. Observe that any solution of (17. 7p has a 
velocity vector of constant length. It is easy to see from the following system 

Uiili = 0UiUiU2 + U1U2U3 + 2UiUsU4 + 3U1U4U5 + AUiU^Uq + . . . , 

U2U2 = —IU1U2U3 + OU2U2U4 + IU2U3U5 + 2U2U4UQ + . . . , 
(7.8) Ms-its = —2uiUsU4 — IU2USU5 + OU3U3UQ + . . . , 
U4U4 = —3U1U4U5 — 2u2U4Ue + . . . , 
U5U5 = -Auiu^ue + . . . , 
UqUq = ... 

Then, 

71 12 " 

a\u\ - 

— — = 2 ^{ukUk + UkUk) = 0, 

k=l 
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for any n, thanks to the cut form of our vector fields and the skew symmetry of fl7.8p . 
The simplest solution may be deduced for constant driving terms Uk, k = 1, . . . ,n. 
The Hamiltonian system immediately gives the geodesic 

Ci = Mi(0)s + Ci(0), 

C2 = M?(0)s2 + U2(0)s + C2(0), 

3 2 2 

C3 = 3Mi(0)(M?(0)y +M2(0)y + C2(0)) +2m2(0)(mi(0)^ + Ci(0)s) +M3(0)S + C2(0), 



In general, Cn becomes a polynomial of order n with coefficients that depend on the 
initial data c(0) and on the initial velocities 'u(O). 

The Lagrangian C corresponding to the Hamiltonian function H can be defined 
by the Legendre transform as 

n k—l ^ n 

c = {c,ij) - H = (Jki^k + i^k ^{j + i)c/A.-i) ^ 2 ^ '^'^■'^^ 

k=l j=l k=l 

Taking into account that 

k-l 

4'kCk = + i)cj4}Jk-j + i'kk- 

Summing up over fc, we obtain (c, ■?/') = Yl^=ikk = Ylt=i'^kUk, that gives us 

1 " 

k=l 

All these considerations can be generalized for n — oo. Thus, we conclude that 
the coefficients of the function p{z, t) in the Lowner-Kufarev PDE play the role of 
generalized moments for the dynamics in Ain and Ai with respect to the Kirillov 
basis on the tangent bundle. Moreover, the L^-norm of the function p on the circle 
is the energy of such motion. 

8. SLE and CFT 

In this section we review the connections between conformal field theory (CFT) 
and Schramm- Lowner evolution (SLE) following, e.g., [B], [13]). SLE (being, e.g., 
a continuous limit of CFT's archetypical Ising model at its critical point) gives an 
approach to CFT which emphasizes CFT's roots in statistical physics. 

SLE^ is a x-parameter family of covariant processes describing the evolution of 
random sets called the SLE;^ hulls. For different values of x these sets can be either 
a simple fractal curve x G [0, 4], or a self-touching curve x G (4, 8), or a space filling 
Peano curve x > 8. At this step we deal with the chordal version of SLE. The 
complement to a SLE^ hull in the upper half-plane H is a simply connected domain 
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that is mapped conformally onto H by a holomorphic function g{z,t) satisfying the 
equation 

do 2 , , 

dt g[z,t)-e.t 

where = y/^Bt, and Bt is a normahzed Brownian motion with the diffusion constant 

2t 

X. The function g{z,t) is expanded as g{z,t) = z -\ \- . . . . The equation (18. ip is 

called the Schramm-Lowner equation and was studied first in |28j [30]. see also |43] 
for basic properties of SLE. Special values of x correspond to interesting special cases 
of SLE, for example x = 2 corresponds to the loop-erasing random walk and the 
uniform spanning tree, x = 4 corresponds to the harmonic explorer and the Gaussian 
free field. Observe, that the equation (18.11) is not a stochastic differential equation 
(SDE). To rewrite it in a stochastic way (following [6], |13j ) let us set a function 
kt{z) = g{z,t) — C,t, where kt{z) satisfies already the SDE 

2 

dkt{z) = -r^dt - d^f 
kt{z) 

For a function F(z) defined in the upper half-plane one can derive the Ito differential 
(8.2) dF{kt) = -ditL.iF{kt) + dt{^L\ - 2L^2)F{kt), 

with the operators L_i = — ^ and L„2 = These operators are the first two 

Virasoro generators in the 'negative' part of the Witt algebra spanned by the operators 
—z^^^-j^ acting on the appropriate representation space. For any state \ip), the state 
L^i\ip) measures the diffusion of under SLE, and (f — 2L_2)\'^) measures the 
drift. The states of interest are drift-less, i.e., the second term in (18. 2p vanishes. Such 
states are annihilated by ^LFii — 2L_2, which is true if we choose the state as 
the highest weight vector in the highest weight representation of the Virasoro algebra 
with the central charge c and the conformal weight h given by 

(6-x)(3k-8) , 6-x 

c = , h = , 

2x ' 2x ' 

and the operators L_i and L_2 are taken in the corresponding representation. It 
was obtained in [6j and |13j . that F{kt) is a martingale if and only if (fi^^i — 
2L_2)F{kt) = 0. We define a GET with a boundary in EI such that the boundary 
condition is changed by a boundary operator. The random curve in HI defined by SLE 
is growing so that it has states of one type to the left and of the other type to the 
right (the simplest way to view this is the lattice Ising model with the states defined 
as spin positions up or down). The mapping g satisfying (18.10 'unzips' the boundary. 
The primary operator that induces the boundary change with the conformal weight 
h is drift-less, and therefore, its expectation value does not change in time under 
the boundary unzipping. Hence all correlators computing with this operator remain 
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invariant. Analogous considerations one may provide for the 'radial' version of SLE 
in the unit disk, slightly modifying the above statements. 

Observe that in this formulation two Virasoro generators can generate the non- 
trivial 'negative' part of the Witt algebra by the commutation relation 

[Lrm — {u — m)L„_|_^. 

There are many forthcoming directions that can follow the described study of 
the Lownere-Kufarev equations. One of possible directions is to study the sub- 
Riemannian geometry naturally related to the bracket generating structure of the 
Virasoro algebra and to the hypoellipticity of the drift operator for SLE. Another is 
to consider analogues of SLE in the case of infinite degrees of freedom (stochastic 
version of the Loewner-Kufarev equation). The alternate Lowner-Kufarev equation, 
infinite dimensional controllable systems analogous to one considered here, are also 
new objects to study. We hope this survey will encourage a new wave of interest to 
this classical subject. 
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